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PREFACE 


This Final Report covers work performed on NASA Grant NSG-1502 between 
1 February 1978 and 31 January 1982 under the Technical Cognizance of 
Mr. Richmond P. Boyden, Transonic Aerodynamics Division, NASA Langley 
Research Center. 

The use of trade names in this paper is essential to a proper under- 
standing of the subject material; their use in no way constitutes 
official endorsement, either expressed or Implied, by the National 
Aeronautics and Space Administration. 
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C3IAPTER 1 
INTRODUCTION 

Magnetic model suspension auid balauice systems for wind tunnel use (1) 
have been designed, tested amd used at MIT's Aerophysics Laboratory for over 
eighteen years. These research programs which have been funded by the USAF, 
NASA, Army Ballistics Research Laboratory and the Army Research Office among 
others have led us to explore the many unique capabilities of the magnetic 
balance for aerodynamic testing. Among these are the precise measurement of 
forces (2,3,4), the measurement of wakes behind cones (5,6), spinning bodies 
at angle of attack (7) , as well as forces and moments on spinning bodies (9) 
cuid studies of ring airfoils (10) , and the production of forced simultaneous 
spinning and coning model motion in the subsonic flow (11) . 

Since the original studies of the roll control problem for symmetric 
magnetic suspension by Stephens (12,13) who described possible methods for 
extending control to the roll degree of freedom, only forced roll rate had 
been produced on a consistent basis (7,9,10,11) before the present research 
grant. This research was undertaiken with the goal of demonstrating closed 
loop control of the roll degree of freedom on the NASA prototype magnetic 
balance at the MIT Aerophysics Laboratory, thus showing feasibility for a 
roll control system for emy large magnetic balance system which might be 
built in the future. During the research under this grant, study was directed 
toward the several aoreas of torque generation, position sensing, model con- 
struction and control system design. These efforts were then integrated to 
produce successful closed loop operation of the analogue roll control system. 
This experience indicated the desirability of microprocessor control for the 
angular degrees of freedom and work was started in this direction. This 
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phase of the work, however, could not be completed within the available 
budget, so work was stopped. Each of the three parts, roll position sensing, 
production of roll torque and stabilization of the feedback loop will be 
discussed in detail below. 
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CHAPTER 2 

GENERATION OF ROLLING MOMEHT 

Production of rolling moment in magnetic balance systems has been in- 
vestigated by Stephens (12,13). He suggested three possible methods: 

a. Interaction between transverse D.C. field components and a 
treuisverse or nonsymmetrical model magnetization. While this method appears 
to have the greatest torque-producing capability, the presence of the full 
matrix of magnetization and applied field conponents with the resulting 
complicated interactions requires a full digital magnetic balance control 
system (14) to provide data reduction and model control with acceptable 

effort and accuracy. Model construction is also difficult for this case. 

/ 

b. Interaction between a single phase transverse A.C. field of 
controlled transverse cuigle with a model supporting a conducting plcuie or 
loop which cam be located inside the wing area. This method has the ad- 
vantage of producing reasonable torque at low roll rates and is independent 
of the lift, slip, and drag degrees of freedom. This method was used in 
initial roll experiments (12) . 

c. Interaction between a rotating two-phame tramsverse A.C. field 
and a copper plated or other conducting model. This method has been used 
for some time (9,15) to produce high roll rates for magnus testing with the 
NASA prototype balamce, and more recently with spinning amd coning models 
(16) . This method produces higher roll rates than method b, but lower 
torques. It will, of course, also work with the loop described in (b) . 

This method is independent of all other degrees of freedom. 

Because of the need for relatively large rolling moments for testing 
winged bodies, amd the need to use existing balance coils, research reported 
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with the roll system concentrated on method b. A qiialitative model of the 
torque generation process is discussed below. 

Torque on Circular Loop 

Consider a circular loop of wire of radius a, as shown in Figure 1. 

Here an A.C. magnetic field of euigulax frequency u and ait^litude inter- 
sects the plane of the loop at an angle 9. The resultcuit field near the loop 
is the sum of the applied field and the induced field B produced by the 
current circulating in the ring. The field of this current has been given 
by Stratton (17) . In the pl^u^e of the wire loop the field is purely axial 
and Stratton's expression reduces to 


H 


Z 


27t(i + g 


2 

KOc) + E(k) 

^-5 


where K and E are the complete elliptic integrals of argument 


a 


A plot of the field distribution is shoxvn in Figure 2, using tabulated values 
of K and E (18) . Here solid points and the solid line show the field of the 
wire loop alone. When this loop ciirrent is generated by an applied uniform 
field/ B^ sin ut/ the current will build up to a level such that the area A 
equals the eurea B/ leaving zero net flux linking the loop if it has zero 
resistance. 

If the loop has finite resistance, area B will be greater thaui area A 
by a sufficient amount to provide EMF to drive the ring current. Since 
/ Eds about the loop » (d/dt ) (/b • dA) over a loop area, the difference 
between areas B euid A (uncancelled driving flux) will decrease as frequency 
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is increased. There vill also be a phase difference between the applied field 
and the induced field. It is these effects as well as the complex geometry 
of actual models which makes detailed calibrations necessary. When the 
applied field makes an angle 0 to the pleme of the ring, only the perpen- 
dicular component B is cancelled by the induced current, where 

B 3 B sin (u sin 9 . 
o 

The force on the wire acting normal to the plane of the loop produces 
rolling moment. This force is proportional to the current I and the parallel 
field component B^ sin ut cos 9. The torque (rolling moment) acting on the 
loop is found by integrating around the wire loop, resulting in 

T = c a^B ^ sin 29(1 - cos 2 cot) (1) 

o 

where c is a constauit that depends on loop geometry resistauice and frequency. 

As was expected the rolling moment peaks at 9 = 45° and has both average and 
periodic con^onents. From measurements of the torque on a copper ellipse 
of 6.35 cm semi-major axis, 2.22 cm semi-minor axis and 0.3 cm thick in the 
400 Hz magnetic balance field, the constant c was found to be approximately 

0.26 dyne-cm/cm2 gauss^ (19). 

The circuit used to apply variable amplitude 400 Hz power (or quadrature 

1200 Hz power) to the pitch and yaw coils is shown in Figures 3 and 4. Here 

the D.C. is isolated from the roll amplifiers by series capacitors and parallel 

L-C traps are used for isolation of the 20 KH position sensing signals and 

z 

400 Hz (1200 Hz) roll power. Methods for developing the control signals for 
the roll power supplies are discussed later in this report. 
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CHAPTER 3 
POSITION SEIISING 

Many methods of remotely sensing, model position withiii a magnetic 
balance have been used. The electromagnetic position sensor (EPS) employed 
with the MIT-NASA balemce has been one of the most successful because it 
presents the following advantages: 

a. Five components of model position have previously been derived 
very accurately over a large r 2 Uige without first order coupling. For 
control purposes the signals are adequate to stabilize a model over the 
entire cross section of the test section. 

b. This method is not sensitive to perturbations within the test 
section such as smoke, pressure probes and balance operators. 

c. The EPS can be made to be insensitive to model geometry. 
Frequently, different models can be suspended without rettming any of the 
control equipment. Changing from ferromagnetic to copper plated models 
requires a simple adjustment of the demodulators. 

Since the ferromagnetic cores used with the magnetic balance until now 
have been symmetric about the wind axis, roll information has not been de- 
tectable by the EPS. Indeed, such symmetry is desiredile. From a functional 
standpoint one should extend the capabilities of the EPS to the measurement 
of roll angles. In this section an experimental method of developing roll 
signals with the existing EPS coils is presented with a description of the 
required circuitry. 

Although the operation of the EPS system has been reported in detail 
by Stephens (13) for the five degree of freedom case, it is instructive 
for the purpose of developing the roll position sensor (RPS) to describe 
this system from a different perspective. 
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The EPS coil arrangement is shown schematically in Figure 5 with a 
sectional view in Figure 6 and the connection schematic in Figure 7 from 
Stephens (13) . The system consists of a pair of Litz-wire wound Helmholtz 
excitation coils euid 14 pick-up coils arranged as shown on the surface of 
a cylinder. The coil parameters are given in Table 1. A detailed 
mathematical analysis of this process is given in Appendix A. The EPS 
is centered at the magnetic center of the balamce. The excitation coils 
are driven at 20 KHz by a high stability audio-oscillator and power-amplifier 
combination. They produce a highly uniform eucial field within the test 
section. The pick-up coils consist of 30 turns of Mo. 30 teflon-coated 
wire arranged to be orthogonal to the desired component of the excitation 
field. Each coil is connected to the one directly across the test section 
and in the saune sense; for example, a D.C. current applied to the H2A-H2B 
combination in Figxire S would produce an additive field along the y axis. 

The notable exception to this rule is the drag coil pair marked AX in 
Figxire 5. This pair is coaxial with the Helmholtz coils and connected in 
opposition. Since the coils are arranged to link no net flux without a 
model pertxirbation, no voltage should be developed in the EPS when a model 
is not present. In practice, small residual couplings between coils exist. 
These cure nulled in the demodulation circuitry. 

When a ferromagnetic model or a diamagnetic model is placed within a 
uniform oscillating field, the field is perturbed (Appendix A) . A diamagnetic 
model can be constructed by copper plating a ferrous model to a depth greater 
than the effective skin depth of the conductor at the frequency of interest. 

The applied field then "sees" only the copper, because currents flow in the 

* Litzendraht conductor, which is called "Litz-wire," is a cable of many 
fine insulated wires that is twisted in its manufacture such that the 
resultant cable has low resistance to both direct and to alternating current. 
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Table 1 


NASA Prototype Balance 
EPS Coil Parameters 


EPS Coil 
Circuit 

DC Resistance^ 
ohms 

Inductance^ 
millihenry s 

Self-capacitance 

microfarads 

Axial 

14.76 

.919 

.00141 

Lateral I 

7.57 

.468 

.00156 

Lateral II 

10.41 

.517 

.00143 

Lateral III 

7.18 

.416 

.00126 

Vertical I 

12.42 

.509 

.00144 

Vertical II 

9.49 

.502 

.00146 

Vertical III 

12.47 

.488 

.00152 


(1) Measvired with Fluke digital VOM. 

(2) Calculated from resonance frequency 

with . 25 yf mica transmitting capacitor. 

(3) Calculated from self-resonant frequency. 
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conductor, which exclude the changing field. Thus the effect of a dia- 
magnetic model is to cancel a portion of the applied field. 

The practical effect of the copper plate is to remove coupling which 
would otherwise exist between the DC support fields and AC position sensing 
fields. When a ferromagnetic ellipsoid is placed within a uniform field, 
the magnetization within the volume is constamt and adds to the applied 
field. If the ellipsoid is diaunagnetic, y is negative and the magnetization 
opposes the applied field. The resultant field is the sum of the dipole 
field due to magnetization of the model and the applied field. For the 
purpose of discussion it will be assumed that the field builds up in phase 
with the applied field for both the ferromagnetic and diamagnetic cases, 
although in practice the magnetization or induced field will be slightly 
phase shifted due to ohmic losses in the model. The field due to the model 
is also assumed to be tied to the model and to rotate with it for small 
angles. It is this effect which is exploited to detect roll position. 

For the axisymmetric case the process of detection has been described 
by Stephens (13). It is desirable to convert the AC voltage at the pick-up 
coils to a DC voltage proportional to model position. This can be done 
ideally by multiplying the pick-up coil voltage by a reference in phase 
with it amd then averaging over a period. Consider now just a pair of 
field lines iss\iing from the model at a particular time t. The applied 
field and the model are shown in Figure 3. If the model is centered, the 
voltage induced in all the coils is zero. Here the model is moved in drag, 
lift amd pitch. The letters next to the lines of flux indicate the voltage 
developed in each coil due to the component of flux perpendicular to their 
surface. The convention here is that is expanding and the area vector 
is positive radially outward. Only the VI, V2, V3 and drag coils aire shown 
here; the results apply equally well to yaw and slip by symmetry. The 
total voltages developed in each coil are given for each configuration. 
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From these diagrams it follows that to first order pure drag signals axe 
measured by AX without lift or pitch coupling. Lift cam be measured from 
V1-V3 and pitch is given by V1+V2 + V3. It is interesting to note that 
V2 contributes nothing to linear model motions. A general statement about 
position sensing in this mamner is that lineau: motions are detected by 
subtracting coil outputs and rotations au:e detected by summing coil out- 
puts. Similair drawings cam be made to show that the coils are only sensitive 
to motions perpendicular to their surface. 

Roll Loop 

The roll position can only be detected by the EPS coils if the model 
has some magnetic asymmetry about its equatorial plame. To this end an 
elliptical copper loop has been fastened to it at some small amgle of 
incidence, as in Figure 9. If a perfectly conducting loop is considered 
by itself amd a field is generated in it to exactly counteract the component 
of perpendicular to it, the RMS radial component of field is then 

B 

B^ = ^ sin 2 9 (2) 

The end view in Figiore 10 shows the effect of the loop alone on the pick-up 
coils as the roll amgle changes. Since all the tramsverse flux leaving 
the EPS reenters it, voltages are developed in each coil due to each flux 
line as described before. The pitch and yaw signals from the V and H coils 
now contain a roll component. These components are found experimentally to 
be proportional to the sine amd cosine of roll angle <p. 

tVhen the loop is placed around a ferromagnetic model at a small angle 
of incidence, several changes take place in the far field measured by the 
EPS. A relatively strong AC transverse field now exists which magnetizes 
the model along its short axis. While the field due to the model alone 



11 


rotates with it in roll, equation (1) implies that the transverse field due 
to the loop vauries as twice the pitch euigle, thus producing coupling between 
pitch, yaw and roll. Since the model perturbation field has a different 
phase for di^u^agnetic and ferromagnetic materials, the pitch and yaw coil 
signals contain a quadrature component which depends on roll angle. These 
can be demodulated separately to give signals proportional to pitch, yaw' 
and roll. In the prototype roll position sensor RPS, the nature of the 
coupling was asstnned to be linear cind attempts were made to null this effect 
by crossfeeding, pitch, yaw auid roll. This method was only effective for 
small roll angles (about ±30 degrees about the set point) . 

Position Signal Demodulation 

The signals received at the pick-up coils cure sinusoidal auid amplitude 
modulated. Stephens (13) has described a method of demodulation and the 
associated electronics. In order to obtedn roll information as well as 
true pitch and yaw, a roll position sensor was devised. It takes the AC 
pitch and yaw signals from the five degree of freedom system after they are 
combined and demodulates them with respect to separate references. Figure 11 
is a block diagram of the roll position sensor (RPS) . A suitably phased 
reference is formed by adding in phcise and quadrat\ire signals. The incoming 
pitch or yaw signals are multiplied by each reference and then low pass 
filtered to obtain an average value. The output stage contains a mixing 
circuit which caui be used to decouple pitch and yaw from roll. Two identical 
channels are required — one for pitch and sin 9 and one for yaw and cos 9. 

Roll Position Sensor Circuit 

The schematic for the RPS is given in Figure 12. Although the function 
of the circuit is similar to the five degree of freedom EPS system, the 
components have little resemblance. Considerable effort has been devoted 
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to reducing the temperature drift and coupling which have caused small 
position errors in the past. The following areas for improvement were 
identified: 


a. Operational amplifiers in the EPS system have insufficient 
gains at 20 KHz. — Open loop gain is described by Jxing (21) to be the 
single most important parameter in describing the ultimate accuracy of an 
operational amplifier. Total hairmonic distortion, output resistance and 
sensitivity to gain variations are each reduced by the factor 


K 


1 

1 + 6a 


( 3 ) 


where A is the open loop gain, and B is the voltage attenuation of the 
feedback network. These concepts are best illustrated by am exaunple. 

The Nexus operational as^lifiers used in the EPS demodulator have an open 
loop gain of 100 at 20 KHz. If a unity gain inverting stage is chosen, 

0 = 0.5; therefore, K = 0.02. The open loop output impedauice is 5000 ohms 
for this amplifier and is transformed to 100 ohms closed loop. This is 
a large output impedance for an operational as^lifier and successive 
stages will load it unless their input impedances are very high. 

b. The slewing rate of the EPS operational amplifiers is not high 
enough to follow 20 KHz signals of 10 volt amplitude .— -The total output range 
cannot be utilized and the signal-to-noise ratio is reduced. The slewing 
rate, SR, required to reproduce a signal of amplitude a is given by 

SR = ao) (4) 

For the case of 20 KHz, 10 volt signals SR must exceed 1.26 V/yS. The 
Nexus amplifiers in the EPS system only offer a slew rate of 1.1 V/yS. 
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c. The phase sensitive demodulators used in the EPS system operate 
only over a small input voltage range before the transformers satiirate. — 
Furthermore, distortion due to the tum-on transient is observed when the 
signal amplitude is small. 

An effort has been made in the design of the roll position sensor system 
to eliminate some of the above deficiencies. The high frequency amplifiers 
are Analog Devices A0518 operational amplifiers. They have open loop gains 
of 65 db at 20 KHz. For comparison purposes, K= 0.001 for the unity gain 
inverter. This yields a factor of 20 improvement in gain performance. The 
70 V/yS slewing rate of this amplifier allows full utilization of its output 
reUige. These amplifiers are very sensitive to capacitive loads, however, 
euid must be suitably isolated. The ring demodulators have been replaced 
with laser trimmed precision monolithic multipliers. They feature less 
than 1 percent total error and better linearity tham the phase sensitive 
ring demodulators. In addition, they are capadsle of operation over the 
full 20-volt rainge. 

The low pass filter is designed to attenuate the high frequency com- 
ponents of the demodulator output aind give a measure of the DC value of 
the signal. It consists of three stages. The first is a passive low pass 
filter with two poles at 1000 Hz. The second adds two poles at 100 Hz amd 
the third has a pole at 48 Hz. The DC gain of the total filter is approxi- 
mately 16 and its attenuation of components above 400 Hz is greater than 
44 db. The design of this filter has not yet been optimized. Possible 
modifications are suggested later. 

The mixing circuit on the output gives a method of nulling the coupling 
between roll and pitch auid yaw. It is possible to add up to the full roll 
output or its inverse to pitch aind yaw and vice versa. 




14 


To test the output of the roll sensor the model. Figure 9, was 
supported in a pivoted Lucite cage and rolled about the long cucis of the 
ellipsoid while the outputs of horizontal roll and vertical roll were 
displayed on an x-y plotter (19) . The resulting maps are shown in 
Figure 13 with pitch auigle of the ellipsoid axis, 0, as a parameter. 

Note that at 0=0 the map is very close to the ideal circle with usable 
control regions in all quadreuits. As 6 is increased, however, the 
symmetry of the output disappears with a region useible for analog control 
limited to about <^ = 30® to <{> = 150®, for 9 = -5® and -10® and <j) = 210® to 
$ = 330®, for 9 = +5® and +10®. This is the effect which limited the 
operating rauige of the analog roll system euid suggested development of a 
digital control system using microprocessors. 
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CHAPTER 4 

OPERATION OF THE ROLL SYSTEM 


From Eq. (1) it is apparent that a time average rolling moment on the 
inclined loop can be generated when the applied AC field makes a non-zero 
angle with the plane of the loop shown in Figure 14. Two possible methods 
of controlling the total field are; 

a. The total field magnitude can be maintzdned constauit and 
rotated as necess6ury by the control system. The average torque is found 
by resolving B into its components, as in Figure 14. 


T 

a 



sin 2 9 


(5) 


where 6 is the angle between the field vector and the loop plane. The torque 


constant K is defined by comparison with the average of equation (1) . 

3tC 

when the field is aligned with the plcine of the loop, there is no rolling 

« 

moment since it links no flux. 


b. The total field amplitude can be modulated while maintaining 
the resultant angle at i45 degrees to the loop plane. Maximum torque is 
therefore available for any given applied field strength. This method 
was tried first because there is a relatively straightforward meeins of 
developing the roll drive signals for a full rotation of the model. This 
method makes calculating torque from cxirrent simpler since 0 need not be 
determined at each data point. 

This method was implemented by Way (19) ; however, he foiind the system 
was unstable to limit cycle oscillations. This was a result of the fact 
that there is no restoring spring constant at zero applied torc[ue. Method a, 
on the other hand, has finite stiffness in roll at zero angulau: error (9 = 0) 
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since the AC field is fully energized at all times, not just when torque 
is being called for by the control system. 

This system was tested by converting the roll drive circuit of Way 
(19) to one which produced the needed constauit amplitude, single phase 
field which could be rotated through a 45 degree range by a compensated 
error signal. For this method a single phase, constant amplitude field B 
is applied trsuisverse to the wind axis at a small emgle 6 to the plane of 
the model loop. The nominal euigular position of the model is (J), as shown 
in Figure 14. 

The angle 5 is generated by the compensator and is represented by the 
compensaced error signal. Referring to Figure 15, 

<{>^=»5+4),jBl=b (6) 


which has components from the inner and outer saddle coils 


B. =» — (cos 0- + sin (})^ ) 


(7) 


Bq =• — (sin - cos 


( 8 ) 


B. =— (cos 5 cos (j) - sin 5sin<fr + sin6 cos (^ + cos6 sin(|» ) 


(9) 


B = — (sin 6 cos 4) + cos 5 sin (p - cos 5 cos (p + sin6 sin 4> ) 


(10) 


For small angles 5 <<1 radian 

B = ^ (a - 66) 

B- = - (6 -H 6a) 


( 11 ) 

(12) 


where 


a = cos (j) + sin 4>, 8 = sin <j> ~ cos 4> 


(13) 



17 


The torque on the loop can be estimated using Eq. (1) to be 

T « |b 1 ^ sin 2 5 (14) 

2 2 

Normally |b 1 = b but since the 5 = sin 5 approximation has been made, 

I B I ^ is given by 

2 2 2 

B- B. + B (15) 

1 o 

2 2 

- ^ (a^ + 6^6^ - 2a66) + ^ (B^ + bV + 2ct6B) ( 16 ) 

Therefore , 

T <= 2b^(l + 6^)6 (17) 

amd the stiffness with respect to 5 is given by 

« 2b^(l + 3 6^) (18) 

Thus the stiffness at the set point is always non-zero and increases 
quadratically for S 0. Furthermore, it is independent of 

for 5 < < 1, « 2b^ (19) 

If b is constant, the system behaves as a lineau: system for small per- 
turbations and therefore cannot limit cycle. As the model moves to higher 
amplitudes. The experimentally observed limit cycle requires the existence 
of a product of a quadratic signal with the torque rate signal. The present 
model is too crude to show this phenomena. 

A photograph of the model used for system development testing is shown 
in Figure 16. This is am ellipsoidal iron core with a copper loop attached 
to it so the loop axis can be inclined slightly to the ellipsoid axis. This 
method of construction was used successfully for an F-16 aircraft model with 
the loop fomuLng the wing leading auid trailing edges and the offset core 


located in the fuselage 
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The circuit used to develop the modulation for the two roll power 
amplifiers is shown in Figure 17. Here the two position sensor outputs 
cure available for selection depending on the model position. The selected 
position signal is compared to a reference and the error signal is fed out 
to the compensator circuit (not shown) . The compensated error signal is 
then supplied in the appropriate sense to the x axes of the two multipliers 
used as modulators. The Y axes are supplied by a stable oscillator auid a 
phase adjusting network used to compensate for the difference in phase lag 
between the two power circuits. 

Compensator 

Compensator circuits essentially similar to the other degrees of 

freedom were used with the addition of 400 Hz notch filters to remove 

pickup. The moment of inertia and torsional stiffness of the open loop 

system was determined by Ramohalli (21) . He found that for this model 

-5 2 

(Figure 16) the moment of inertia J=3.1 x 10 kgm . With the maximum 
•current of 18 amperes T = .00764 N-m from the inner saddle coils alone. 
Since T al^, linecurizing about zero, we have 3 t/ 3I = Hence, is 

determined to be 0.00764/36 = 0.00042 H-m/amp. k,g and 1/T vary, depending 
on which coil system is being considered. 


k^ = 4.07 


k = 4.31 
o 


A control loop for one coil is shown in Figure 18 and the combined system 
for both coils is shown in Figure 19. 

If a series compensator design is used, we have for the inner saddle 
(the outer saddle is similar except for a few parameters) 

K(s) (s + 140) 


KG G H = -r 
c p 2 


s (s + 296.5) (s + 50) (s + 140) 
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A root loctis for this is shown in Figxire 20 emd the series compensator' 
poles are listed below. 


TABLE 2 Inner Coil Loop Parameters 


K 

k 

Poles in 
minor loop 
configuration 

Poles in series 

compensator 

configuration 

(300) 

4.07 

-21.2 + 21. 2j 

-24 +21j 

30^(300) 

2.04 

-14.5 + 20. 5 j 

-35.5, -12 

(30^) (300) 

8.14 

-26.5 + j20 

-21 +37j 


Notice that in the minor loop feedback. Figure 19, the poles eire not 
very sensitive to gain variations. However, we need rather large gains in 
the loops. In the series compensator we do not need the large gains, but 
for gcd.n variations the poles are very sensitive and even go unstable. 

The minor loop has the advantage of being stable at all gains. It was the 
configuration that was successfully implemented. 

System Operation 

The' error signal (after compensation) is fed directly to the summing 
point on the outer saddle chainnel. The error signal is also inverted and 
fed to the summing point on the inner saddle channel. The operating point 
(euigle) is initially set using the position control pots and P^, 

Figure 17. 

The model position reference and response to various excitation phases 
is given below in Table 3. The position sensor outputs as a function of 
model angle are plotted in Fig\ire 21. These are the x emd y components of 
Figure 13 for zero pitch angle. Note that a different output aind sign 
combination is needed for different zero position euigles as indicated by 
the solid and dotted braces in Figure 21. 
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TABLE 3 Model Response to Different Roll Inputs 





Model Position Reference 

0® Model loop horizontal above pivot in front cuid below pivot in rear. 
Increasing angle -*■ Rolling right wing up. 




Model 

Angle 

Inner 

saddle driven at lOV 

135® 

315® 

outer 

saddle driven at lOV 

45® 

225® 

Inner 

+ Outer in phase 

0® 

180® 

Inner 

Outer out of phase 

0® 

180® 
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The system was operated successfully with both the ellipsoidal model 
and the F-16 aircraft model. The reuige of roll angles was about 45 degrees 
for one setting of the control system as explained above. To extend this 
range, efforts were directed to implementation of digital control for the 
three angular degrees of freedom. This is discussed in the next chapter. 
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CHAPTER 5 

INITIAL EXPERIMENTS WITH A MICROCOMPUTER AND THE MAGNETIC BALANCE 

Initial experience operating the magnetic balance system with the 
cinalog roll control loop indicated that initial recommendations of Luh et al 
(14) and Way (19) that digital control for the system be implemented were 
correct. 

Plans were then made to utilize an existing DEC 11-23 microcomputer 
with two additional Motorola 6800 microprocessors for data acquisition 
and digital control experiments . The work was planned in four steps : 

1. Connect existing analog position emd current outputs and position 
set inputs to the MINC 11-23 to provide digital data readout and position 
control on the cinalog system. 

2. Check out system operation and verify choice of components 
while taking data. 

3. Design microprocessor architecture and software for control 
of roll, pitch, and yaw. 

4. Install and test digital roll control. 

Because step (1) required much more effort than anticipated, it was 
only possible to complete steps (1) and (2) of the work. 

The planned system is shown in Figure 22. Here the MINC is used to 
provide long-term monitoring of model position and to quickly read current 
and position data emd then tremsfer it to the disc. Data reduction was 
performed with the balance system off line. To alleviate the need for am 
extremely large fast microprocessor the separate dedicated MC 6800 units 
were selected with one for each degree of freedom. It later appeared that 
all four might not be required so only two were purchased. However, 
termination of the contract did not allow us to install and test them. 
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In step (1) the position aind current outputs of the euialog control system 
for the magnetic balance were connected with two 27-pair shielded cables to 
the Data-Translation A-D Boards at the MINC output. A block diagrzun for this 
system is shown in Figure 23. Here the computer receives inputs of model 
position from the electric position sensor and magnet current from the shunts 
in the coil circuits for lift, drag, slip, pitch and yaw. A series of ex- 
periments were then carried out measuring the pitch damping of the ogive 
cylinder model previously built for spinning and coning tests (16) using 
white noise excitation of the model in pitch (22) . Also some measurements 
of the asymmetric side force on this model when held statically at auigles 
of attack up to 30 degrees were made. 

These results are of interest because they indicate the feasibility 
of measuring the dynamic derivatives without inducing large amplitude 
periodic model motions (which are difficult to induce with power limited 
equipment) using short streams of data and in a noisy environment. 

In the case of the asymmetric side force at high angles of attack 
the same fast Fourier transform softwcire was used to obtain the spectrum 
of the measured side force. It is hoped this data will shed some light on 
the nature of the flow causing this side force. 

Pitch Damping Measurement Using Raindom Excitation 

The procedure followed here was to adjust the model suspension to be 
very stiff in all degrees of freedom except pitch. The pitch position 
response to white noise excitation was then stored during a run and analyzed 
using the MINC. For this case all displacements except pitch are negligible 
and the equation of motion is approximated as: 
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where J = Moment of inertia about the center of rotation of the model 
0.. = Damping due to the magnetic suspension feedback system 

II 

a Stiffness due to the magnetic suspension feedback system 
P(t) =■ Pitch input driving function 

and the aerodynamic forces have been approximated by Mq 9 + M . Note that 

• • # • • 

for this pure pitching motion 9 = a so the term m 9 = M*xa + Mix9. Also, 

C( w 

• • 

since q ( nondimens ional pitch rate) =» iref 9/2U , C„ xq = C„,x9. Thus the 

q 9 

• • ^ 

term M is a combination of a emd 9 terms » Using LaPlace transform we can 

tedce advantage of the fact that the transform of P(t) , P(u) = constzmt if 

P(t) is white noise. In the experiment white noise was obtained from a 

generator connected to the control system. Hence, in the transform domain, 

2 

9 (oj) /const = 1/ TJu + (D -M)o) + (K -Mq) ] . The beuidwidth of the corres- 

M MO 

ponding power spectrum Aa» = 2(u where u = the natural frequency, Cem be 

n n 

used to determine the aerodynamic damping in pitch. 


Ab) 




With the wind off, M * 0 auid we cam measvire D^J. 

J, the moment of inertia, is easily measured either directly by torsional 
pendulian or by measuring the resonant frequency and torsional magnetic stiff- 
ness wind off. 

Thus M Ccui be determined by subtraction. (Note that Mq; and hence C„ 

u M 

can be found in similar fashion from the natural frequency. However, 
sufficient data has not been obtained to determine the accxiracy compared to 


the classical static approach.) 

_ * i M 9 

Then C,, xa + x9 = = 

M» MS 1 2 

a 9 ^ puj SD 


or since a = 9 


C + C 
M* M 
a 


2 ' — 00 
M 


■9 i SD 
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This method of determining pitch dcut^ing is not new; in fact, it is 

the basis of the free oscillation technique used for years in wind tunnels. 

What is new is that the combination of the magnetic suspension system with 

the minicomputer permits rapid acquisition of data with high resolution 

combined with very small random model motion. 

It is also possible to exploit the capability of producing a wide 

variety of model motions. For example, by applying the above technique in 

lift auid side slip, C , C eUid C and C can be determined. This caui 

^8 ^8 

of course be done about zero degrees or about any desired angles of pitch 
or yaw. 

Pitch Damping Data 

The body tested was 2 ui ogive cylinder of 1 inch diameter and 5 inch 
length with a 1.506 caliber tcingent ogive nose of 2.519 radius. This was 
the same body for which Magnus measurements were reported by Birtwell (9) . 
FFT's of the wind off, wind on and digitally filtered wind on data are 
shown in Figures 24, 25 and 26. These were processed off line from data 
records containing 256 points of data recorded in one second. System 
response ctirves resulting from applying Ramohalli's data reduction approach 
(22) are shown in Figures 27 and 28. The resulting values for the dan^Jing 
coefficient (C. + C.. ) aure given in Table 4 as a function of wind speed for 

a white noise amplitude equivalent to 0.2 degrees pitch angle. 

Data was taken on two model cores. One was a soft iron core amd the 
other was a permamently magnetized samairiuitt-cobalt core. This method was 
used to determine if the measured damping was being affected by hysteresis 
in the iron model. As cam be seen from the data the model with the two 
cores gave results which were quite different. These results are shown in 
Figures 29-33 for the iron core and Figures 34-40 for the samaurium-cobalt 
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permanently magnetized core. Figure 41 is a coit^>arison cross plot of the 
200 fps data versus angle of attack. These results differ so much that a 
shift in the control system induced magnetic deunping from hoiir to hour is 
suspected. The most likely cause of this is inconsistent location of the 
analog system poles and zeros resulting from noisy connections or drifting 
con^onents . 

Because this method of determining dynamic stadsility derivatives 

depends on subtracting out the wind off response, any change in system 

performance function between wind off and wind on conditions is reflected 

directly as an error in the measurements. This initial experience indicates 

that the addition of a computer to the magnetic balance has greatly improved 

system versatility euid data rate. The digital computer mcdces possible use 

of white noise excitation at small amplitudes for determining dancing in 

pitch with much shorter run times than preceding methods. The technique 

should also be applicable for measurement of C,. and C, . from a pure 

a a 

plunging motion and c emd c from a slipping oscillation. Since the 
power demands cure much lower than for forced oscillation testing, system 
operating range should be wider for this method of testing. 

It has also demonstrated a need for greater stability of the analog 
control circuits before dynamic data accuracy can approach the level 
attained in static measvurements . Digital control may be the best approach 
to improve stability. More effort is needed to develop this technique and 
also to apply the digital data system to determination of dynamic stability 
derivatives by direct instantaneous force measiirement rather than by 
pcurameter identification methods. 

While data was being taken at high angles of attack, some static 
measurements were made of the side force at 30 degrees angle of attack. 

These results eu:e presented in Figxire 42. 
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TABLE 4 Prelimineury Pitch Damping Data for Ogive-Cylinder 
Model 


Riin 

Angle of Attack 
a (degrees) 

V 

ft/sec 

Soft Iron Model 

2?0J 

n 

C + C 
M M 

q < 

/deg 

Tcure 

2 

0 

66.5 

— 

1 

2 

50 

71.4 

-665 

2 

2 

100 

73.7 

-489 

3 

2 

200 

69.5 

-103 

4 

2 

300 

71.5 

-113 

Teire 

4 

0 

70.9 

— 

1 

4 

SO 

65.6 

t713 

2 

4 

100 

71.9 

-67 

3 

4 

200 

86.6 

-532 

4 

4 

300 

67.8 

+71.6 

Taure 

8 

0 

80.0 


1 

8 

50 

77.1 

+393 

2 

8 

100 

70.2 

+666 

3 

8 

200 

64.9 

+514 

4 

8 

300 

74.3 

+131 

Tare 

16 

0 

73.6 

— 

1 

16 

50 

66.2 

+999 

2 

16 

100 

64.2 

+637 

3 

16 

200 

68.1 

+186 

4 

16 

300 

68.7 

+111 

Tare 

32 

0 

73.4 

— 

1 

32 

50 

83.2 

-13C0 

2 

32 

100 

79.7 

-433 

3 

32 

200 

99.6 

-893 


Tare 

2 

Samarium-Cobalt Model 
0 

79.0 


1 

2 

50 

98.7 

-2670 

2 

2 

100 

73.2 

+392 

3 

2 

200 

66.3 

+432 

4 

2 

300 

150 

-1610 

Taire 

4 

0 

120 

— 

1 

4 

50 

111 

-1260 

2 

4 

100 

121 

-54.4 

3 

4 

200 

119 

+42.2 

4 

4 

300 

58.2 

+1411 
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TABLE 4 (continued} 


Run 

Angle of Attack 

V 

2 SO) 

n 

C + C 

M' 


a degrees 

ft/sec 


/deg 

Tare 

8 

0 

105 

— 

1 

8 

50 

95.4, 

-1360 

2 

8 

100 

93.0 

+842 

3 

8 

200 

133 

-936 

4 

8 

300 

98.9 

+146 

Tare 

16 

0 

109 

— 

1 

16 

50 

73.0 

+4820 

2 

16 

100 

87.2 

+1450 

3 

16 

200 

101 

+250 

4 

16 

300 

104 

+113 

Tare 

20 

0 

107 

— 

1 

20 

50 

83.5 

-3150 

2 

20 

100 

83.3 

+1590 

3 

20 

200 

62.6 

+1330 

4 

20 

300 

71.7 

+793 

T 2 u:e 

24 

0 

77.4 

— — — — 

1 

24 

50 

109 

-4280 

2 

24 

100 

76.3 

+69.3 

3 

24 

200 

117 • 

-1340 

4 

24 

300 

73.5 

+87.9 

Tare 

28 

0 

84.3 



1 

28 

50 

95.4 

-1500 

2 

28 

100 

71.9 

+847 

3 

28 

200 

111 

-913 


Note: No transition devices were used on the model 
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Chapter 6 

Summary and Conclusions 

During the research under this grant, study was directed toward the 
several areas o£ torque generation, position sensing, model construction 
and control system design. These efforts were then integrated to produce 
successful closed loop operation of the analogue roll control system. 

This experience indicated the desirability of microprocessor control for 
the angular degrees of freedom and work was started in this direction. 

It progressed through: 

1. Connecting existing analog position and current outputs and 
position set inputs to the MINC 11-23 to provide digital data 
readout and position control on the analog system. 

2. Checking out system operation and verifying choice 
of components while taking data on pitch damping. 

This research indicated that: 

1) The electro-magnetic position sensor (EPS) could be expanded to provide roll 
angle position information suitable for closed loop model control 

x^ich was achieved. 

2) The txjo phase AC method of torque production was capable of 
supporting an F-16 type model at aerodynamic loads consistent 
with the rest of the system. The available roll moment was limited 
by two factors a) high AC resistance in the saddle coils, 

b) limited roll power (2 kw/ channel). 

3) The EPS roll position output was a sufficiently complex function 
of model roll, pitch and yaw angles that microcomputer control is 
desirable for static tests and mandatory for dynamic tests. 
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4) The MING 11-23 micro-computer used for pitch acquisition and model 
position set provided much better data resolution and about 10 times 
faster data rate than the analog system. 

5) Additional work is needed to make control system stability and 
accuracy consistent with this new level of data resolution, parti- 
cularly for dynamic tests such as pitch and yaw damping or measure- 
ment forces on a model undergoing complex motion such as spinning 
and coneing . 

Recommendations 

It is recommended that the original plan for implementing digital 
control in pitch, yaw, and roll be completed and additional computer 
facilities be added to permit the new range of dynamic testing to be 
explored. 

The magnetic suspension and balance system provides the unique 
capability for sting-free testing of models undergoing complex motion. 
This appears to be the best method for measuring the results on flight 
vehicles of the complex aerodynamic phenomena associated with these 
motions. This capability should be exploited. 
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Schematic of E.P.S. transducer coils 


Figure 7 
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Figure 8 


EPS signal contributions from single flux lines 





Figure 9 Radial Component of Field 
from a Conducting Loop 
in a Uniform Field 



Roll Signal Generation 
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Figure 10 


Figure 11 ROLL POSITION SENSOR BLOCK DIAGRAM 

TWO REQUIRED 


IN PHASE 
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Figure 12 


ROLL POSITION SENSOR CIRCUIT 
(ONE OF TWO REQUIRED ) 
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Figure 13 x-y Plotter Display of Roll Sensor Outputs as 
Model is Rolled at Various Pitch Angles 


Figure 14 


Transverse Field applied to 
a Model with a Roll Loop 
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FIGURE 15 resolution OF FIELD COMPONENTS 



Figure 16 Ferromagnetic Model with 

adjustable Copper Roll Loop 
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Figure 18 Block Diagram of Control Loop 
(2 req'd; 1 per coil) 
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Figure 19 Control Loop for Both Coils with Proper Constants 
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Figure 20 


Roll Control Root Locus 
Series Compensation 


Vertical 



Figure 21 Roll Position Sensor 
Output Signals 


Model Angle 


Model Angle 
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Figure 22 Block Diagram of Proposed Microprocessor 
Control System 
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Figure 23 Block Diagram of Computer-Magnetic 
Balance Combination 
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Figure 24 FFT of wind on noise off data. 
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Figure 26 FFT of wind- on noise 
peaks at 60 and 120 
of the filter. 
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Figure 28 FFT of the analytical model of the wind on data 
for O.J.*pitch noise. The broadening of the peak 
is noticeable. (cf Fig. 27) 
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Figure 33 Pitch damping vs velocity for iron model at 
32 degrees incidence 
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Figure 34 Pitch damping vs velocity for samarium cobalt model 
at 2 degrees incidence 





Figure 35 Pitch damping vs velocity for samarium cobalt model 
at 4 degrees incidence 
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Figure 36 Pitch damping vs velocity for samarium cobalt model 
at 8 degrees incidence 





Figure 37 Pitch damping vs velocity for samarium cobalt model 
at 16 degrees incidence 
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Figure 38 Pitch damping vs velocity for samarium cobalt model 
at 20 degrees incidence 
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Figure 39 Pitch damping vs velocity for samarivim cobalt model 
at 24 degrees incidence 




69 



Figure 40 Pitch damping vs velocity for samarium cobalt model 
at 28 degrees incidence 





Figure 41 Pitch damping vs angle of attack for samarium cobalt 
model (hatched) and iron model (open) at 200 fps 




Angle of attack = 30 degrees, angle of yaw = lo degrees. 


Figure 42 Side Force Coefficient vs Wind Speed 
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Appendix A 

ESTIMATION OF MOTION-INDUCED CHANGE IN COUPLING 
BETWEEN A MODEL AND A PARTICULAR EPS COIL 

1. It is well known that the induced EMF in a pick-up coil 
may be deduced from the integral 
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Hence it is sufficient for these purposes to compute the 
normal component of the field as it is effected by the 
model and core geometry. Note ; We assume the 
magnetization is instantaneous; i.e., ^ = f (geometry) 

dT 3t 

However, a lag may exist. This lag is of the order of a 
few seconds, which is very long compared to the time required 
to establish the "magnetizing current". The 20 KHz field may 
thus never reach the shape of the D.C. field. 

2. The calculations will be carried out in the following way; 

a. The magnetic potential for a sphere will be found 

to within an arbitrary multiplying factor that 
represents the uniform 20 KHz exciting field. 

b. The magnetic potential will be expressed in a 

cylindrical coordinate system so the pick-up 
coils consist of four elements: two elements at 

constant radius (r) and constant (x) , different 
x's; and two at constant (r) and constant angle (|J>) . 
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c. The radial component of this field will be com- 

puted and integrated over the area of the coil. 

d. When the model is inclined at some angle 0/ the 

second potential is needed to describe magnetiza- 
tion due to the field in the normal direction; 

i .e . , 


H=HX 


,cos^0, sin^0^ 
i-R — + -R ) 1+ 


(sinOcOsG) (j 


1 


Thus we allow^for the eccentricity of the model 
in computing H even though we have assumed a model 
so small the induced magnetic potential is 
essentially that due to a sphere. A similar state- 
ment can be made for a general ellipse. 

The exact solution for the induced magnetic 
potential by an axially symmetric ellipse is 
known. The expansion of that potential in 
spherical coordinates is given in Appendix B, 

The result is given below in Paragraph 4. 


3. A sphere of radius a centered at the origin, with 
E>ermeability ii*in a xaniform field F , has a magnetic 
potential ^ 



O' - 


y+2 



3 

a X 


2, 3/2 
+ r ) ' 


) 


in the x,r$ coordinate system. If z = r sin 


z 


F (r sin $ - 
z 


U-1 a r sin 4> 


M+2(x^ + r^) 


3/2 


) 


In use the uniform field part must be subtracted from Y for 
the cross potential. Here we define ^ 


★ 

Note, the case of a perfectly conducting sphere (or approxi- 
mately when the sphere is copper plated several skin depths 
thick) can be found by setting y = 0. 



A 




and P 


'u M 


H sinGcosO 
■^C 


M 


z 



Since we are interested in the effects of radial motion/ two 
alternative procedures present themselves. In the first we 
offset the center of the magnetized body from 0/0 to Ar/ 

$ The distance in the denominator of the magnetic 

o • 

potential thus becomes 

+ Ar^ - 2rAr cos ($ - 


and the derivative of this terra/ 


i_( 

3r 


-) = 


2 (r-Ar cos($-$ )) 
o 


X +r +Ar -2rArcos ($-<!> ) 

o 


(x^+r^+ (Ar) ^-2rArcos 


The second procedure is to change the position of the coil 
and is much more complicated. 

If the model were to have a yaw angle/ / the perturbation 
field in yaw becomes 


SAY = M 

y y 


(- 


y-l a r cos <f> 


y+2 (^2 ^ ^2j 3/2 


Here M 

y 


is the magnetization in the yaw direction if it exists. 


Note: S = - V4'. 
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4. In ellipsoidal coordinates u,v,$, one finds 


Y = F 


cosh u cos V + 


|((J) -1) (y-1) cos V 

G 


^ ,1 |_j„ 1^, 

c 


/ •. . T- u cosh u- ^ 
( 1 + cosh 2 to ) 


where a = semimajor axis (in x direction) 
b = s emimino r axis 

c = /a^ - 

a = c cosh u , b = c sinh u 
o o 

X = c cosh u cosv 


r = c sinh u sinv 
<I> = i (tan"^ -) 

y 


If this is expanded into a form compatible with spherical 
coordinates in powers of (c/a) (c -»■ 0 for a sphere; one obtains 


y = F (x) 
x " 


5. The coils of the EPS which are of interest can be 
described by two straight lines and arcs of two circles. 
Thus we are interested in the area integral 


I 



M (x,R,$) dx Rd$ 


Here R is the radius of the cylinder of the coil axis. 
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f 

^r^+(4r) 


6. With all this background we can now compute the radial 
component of the field 


H 


r 


3Y 

3r 


H =H •(3 
r x) y+2 


a\ . 


(r - Ar cos (<5-$^)) 


(x^+r^+ (Ar) ^-2rArcos 


oii \ 

5/2 C 

-$ ) ) 


+ sin 0 cos 0 ^ a^ (1 - ^ 


2 2 2 

Arcos ($-$^) ) sin$-3 (r +(Ar) -2rArcos ) (r-Arcos ) sin $ 

^-2rArcos ) (x^+r^+ (Ar) ^-2rArcos ($-$ ) ) ^^^r (x^+r^-2rArcos ($-$ )) 

o 00 



Note, 


(x^+r^+Ar^2rArcos 


, 2^ 2\n/^ 
(x +r ) 


(1 + 


nrAr 

2^ 

X +r 


„ cos (<!>-$ ) + . . 
2 o 


This latter approximation neglects terms 


0 ( 


(Ar) 

2^ 

X +r 


2 

2 
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7. Making the linearized assumptions 


~ M+ 2 ^x 2 +rV /2 ^cos($-$^)) 


2 

x~+r‘ 


1+ ^ ^ cos 

+ sinOcosOd - — )sin$( y 2 ^ ' I, -3 / 2 

C (x +r ) ' 


3r^ (1-5(1 - ^ cos(»-t„) 


_r — cos(»-»^) 

X +r V 

, 2^ 2,5/2 \ 

[x +r ) ' J 


8. Comparing Paragraphs 5 and 7 we see there are eight 
integrals involved in approximating e in Paragraph 1. 


*2 ='2 


*2 X2 


_ 2 f f 3x dx do . _ „2 r r 

1 -^ J J (^2^r2,5/2' 2 - ■* J J 

*1 ='1 *1 Xi 


SxcosC^-iI* ) dx d$ 

(x 2 +r 2 ° 5/2 


I, = R 


*2 ='2 

_fi 

*1 ='1 


15xcos ($-<!> ) dx d<I> 

(x2+r2°7/2 ' 


I. = R 


^2 ^2 

S( 

x^ 


dx d$ sin $ 
(x2+r2)3/2 


I. = R' 


^2 ^2 

(I 

$1 x^ 


3cos($-$^)dx d$ sin 4> 
(x2+r2, 5/2 '• 


I- = R' 


*2 Rj 

n 

♦1 ='1 


3 dx d$ 


2 2 
(x^+R^) 


sin $ 
5/2 
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= R' 


*2 '‘2 

J/ 

*1 


15 cos($-$^)dx d$ sin $ 


( x ^ R ^)=/2 


“s' 


I„ = R 


j f “ 


15 cos($-$ ) dx d$ sin $ 


9. ThuS/ 


^ sinOcosSd - % (l^-l 5 +( 6 l 5 -l 5 )^)> 


This relation is of great value because it allows us to estimate 
the performance of a given coil in response to motion of the 
magnetized model. 

I. Centered model 



3 

a 


y-1 

y+2 


R‘ 


<Vh> 



+ r2)V2 



+ r2 )V2 


e is the reference value. It may vary from coil to coil, de- 
pending upon how well R, 4>, x are known. Note e- is propor- 
tional to the volume of the model and the area 8 f the coil. 
Note, too, the perturbation field decays like the reciprocal 
of distance cubed, which is quite rapid. 

II. Model perturbed a distance Ax, angle of attack zero 


e = e -M a^ ^ (i> -<I> ) ( 

®Ax o x^ y+2 " '^2 ^1^ ^ 


3Xj^R 





3X2R 


Ax 

R 
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III. Model perturbed in radial direction, Ax=0, angle of 
attack zero. If = (o,tt), this is y displacement, if 

$ =(■?/ -^) / it is vertical 
o z z 


®Ar$ 

o 




M 

R 


Where 93-93 = 3/2 ~ (^^2^r2,3/2 " 5/2 + 


+R ) 


IV. Model pitched an angle 0, Ax=0, Ar=0 


e = e_+M..a^ * R^sin0cos0 (1- (sin $~-sin $T){g4-0c) 


o X 


y+2 


where 94“9g“(l ” 


R‘ 


R^x 


3R 


r) ( 




2 2 






2 2 
+R 


2 2 2 
+ R^)^ 


R^x^yx2^+R ^ 

3(x ^ + R^)^ 


The term in 0 is due to the rotation of the magnetization vector. 
Hence, there is a similar term in angle of yaw when 0=0. This 
constitutes the complete collection of linear terms. 

10. Cross coupling terms abound. Clearly from Paragraph 7 
one can see Ax, Ar, A$ products. There is a suppressed 
coefficient of the 3xr, etc. term, namely 

2 2 
cos 0 , sin 0 


which must be included if second order terms are to be accounted 
for. The formula in Paragraph 7 has one other implication. If 
<I>Q = ti / 2 , the perturbation is in the z direction 
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Ar cos sin 4> 


This moans three coils are necessary to sort out Ax, Az and 
0 to a first order because there are three terms that con- 
tribute to AM^ due to displacement. That is, the system 

tends to solve throe simultaneous equations for each set of 
displacements. Note a redundancy exists in Ax! This re- 
dundancy is removed if roll is measurable because now one 
must solve all the equations at once; i.e., the magnetic 
potential may be written 


a^x 


4* u+2 * 2 2 2 3/2 

X y+z (x^+r +(Ar)^ - 2rArcos ($-<!> )) ' 


M.. ( 


i^v^^+(Ar) ^-2rAr 


cos ($-'1'^) cos 


^ (x^+r^+ (Ar) ^-2rArcos 


- M 


y-1 
z y+2 


( 


a^/r^+ (Ar) ^-2rArcos sin $ 

(x^+r^+ (Ar) ^-2rArcos (‘I’-'I*^) ) 


Now, ^x'^y'^z functions of pitch, yaw and roll angle 

through the usual transformations. This implies three con- 
tributions. The X, Ar, implies three additional contributions. 

Hence, six coils are needed. 

In the past the sensitivity was doubled by adding or subtracting 
opposed coils in easily defined ways. Thus, the number of coils 
becomes six for plane motion and only twelve for complete motion! 
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Appendix B 


Table of Integrals 


*2 ='2 

1 =“' J I 
*1 *1 


3x dx 




/ 


2x dx 




= R ( 


(X 3 ^ 2 ^r 2 )V 2 ,x32h.r2,3/2' 


l2=R 


*2 ="2 

/i 

=‘1 


3xcos($-$ ) dxd$ ~ 

uW)°/ ~^ = R (sin(<.2-<f„)-sln(t^-J^)) 


(x^2+r2)V2 ,^^2^^2,3/2' 


l3=R 


*2 >=2 

// 

*1 X 3 


15xcos($-$ ) dxd<t . 

= y - y% = 3R"*(sin($ -$ )-sin($,-$ )) 

(x^+r"') ^ ° ° 


‘(x, 42 )=/ 2 ‘ (x,W) = / 2 > 


I* = R 


*2 ==2 

// 

=<1 


sin $dx d<I> dx 

(x2-hr2)3/2 


R(sin ‘J> 2 ~sin 


’■/ 


dx 

(x"+r2)V2 
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All these integrals are in the form f(y)*g(x). 
Thus, = Rf^($)*g^(x) 


where f^{<J>) = sin$ 2 ~sin«I>j^ 



2 . 2 
— cos $_+COS 

Note ^5 = ^4 and f^=f^=fg=cos ( 2 sin(j>^( — 


sin$ 


(sin 2 $2 “ 2 


y dx 

— j — w" c From Pierce Integral 

(x +R^) 


Table, page 24, #164 
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/ 


dx 




xl/ x^+R^ 

3 R^(x^+R^)^ 


2 2 

/•*/ 


dx 


(x^+R^) /x^+R^ 


But the second integral is clearly g^. 


Similarly, 


^8 = 


i 


dx 


_ X Jx^+R^ 

(xW)3.^?;^"5 r2(x2+r2) 


A ^ 

I -s-f-r 


dx 


(x2+r2)2^:? 


= 


X .CW 

5R^(x^+R^) ^ 


5R 


2 ^5 


9c: = 




3R^(x^+R^) ^ 


3R 


2 ^4 


= 


„2 i 2. „2 7 ) 2. „2 

R •'^2 ^ *^^1 


^5 ^5^5' 

I, = 15R^ f^g^, 


^6 = 3r 3 £,95 
Ig = 15R fgSg 
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